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Abstract 

The shape of the lepton spectrum in inclusive semileptonic B — ► X £ T>i decay 
is sensitive to matrix elements of the heavy quark effective theory, A and \\. 
From CLEO data we extract A = 0.39±0. 11 GeV and X x = -0.19±0.10 GeV 2 , 
where the uncertainty is the la statistical error only. Systematic uncertainties 
are discussed. These values for A and Ai are used to determine \V c b\ and the 
MS bottom and charm quark masses. We discuss the theoretical uncertainties 
related to order (AQCD/^b) 3 effects and higher orders in the perturbative 
expansion. 
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The operator product expansion (OPE) shows that in the limit m& ^> Aqcd inclusive 
semileptonic B decay rates are equal to the perturbative b quark decay rates [|1|. Experimen- 
tal study of such decays provide measurements of fundamental parameters of the standard 
model, such as the CKM angles \V c b\, \V u b\, and the bottom and charm quark masses. 

To obtain precise theoretical predictions for inclusive semileptonic B decays, it is im- 
portant to be able to compute nonperturbative effects suppressed by powers of AqcdA?V 
There are no nonperturbative corrections at order AqcdA^&j arid the corrections of order 
(AQco/^b) 2 are characterized by only two matrix elements [0|§ 

Ai = (M(v) | K m 2 K I M(v)) , (1) 
2m M 

and 

A 2 = — (M(v) | h v 9 - h v | M(v)) , (2) 

ZdM m M & 

where M is a pseudoscalar or vector meson containing a heavy quark Q, and h v is the 
heavy quark field in the effective theory [Q with velocity v. du = 3, — 1 for pseudoscalar or 
vector mesons, respectively. The decay rates also depend on the quark masses, which can 
be expressed in terms of the heavy meson masses and the parameters Ai, A2 and A, where 

1 t -^1 + ^mA 2 , , 

m M — m Q + A — ^ +..■■ (3) 

The matrix element A 2 is then determined from the measured B* — B mass splitting, 
A 2 — 0.12 GeV 2 . The quantity A || also sets the scale for the deviation of the exclusive 
B — > D(*> ivi decay form-factors from the Isgur-Wise function 0. The analogue of A in 
the baryon sector, Aa = m^ b — m& + . . ., describes all Aqcd/'^q corrections to A& — > A c £ve 
decays 0, and is related to A via Aa = A + — m B + . . .. 

To carry out accurate calculations it is crucial to have reliable determinations of A and Ai. 
In the past, these quantities have been estimated using models of QCD ||, and extracting 
them from experimental data was attempted P|-|i"2"|| . Sum rules were also derived to constrain 



Ai [131, however, perturbative corrections weaken these constraints [ 14| . In this letter we 



extract A and Ai from the shape of the inclusive B — > Xlv^ lepton spectrum, and also 
translate our results into a determination of |V^|, and the MS masses m^mj), and m c (m c ). 

The CLEO Collaboration has measured the inclusive B — > X £ lepton spectrum both 
by demanding only one charged lepton tag ]TJ], and using a double tagged data sample fl6 



where the charge of a high momentum lepton determines whether the other lepton in the 
event comes directly from semileptonic B decay (primary) or from the semileptonic decay 
of a B decay product charmed hadron (secondary). The single tagged data sample has 
smaller statistical errors, but it is significantly contaminated by secondary leptons below 
about 1.5 GeV. For our analysis we use the data as tabulated in Ref. fl7|l . 

The OPE for the lepton spectrum in semileptonic B decay does not reproduce the physi- 
cal lepton spectrum point-by-point near maximal lepton energy. Near the endpoint, compar- 
ison with experimental data can only be made after sufficient smearing, or after integrating 
over a large enough region. The minimal size of this region was estimated to be around 
300 — 500 MeV || . This, and the fact that the experimental measurement of the lepton 
spectrum is precise and model independent only above about 1.5 GeV, impose a limitation 
on what quantities can be reliably predicted and compared with data. On the one hand, we 
want to find observables sensitive to A and Ai; on the other hand, we want the deviations 
from the b quark decay prediction to be small, so that the contributions from even higher 
dimension operators in the OPE are not too important. The observables we use should not 
depend on \ V c b\. Thus we consider 
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Before comparing the experimental data with the theoretical predictions for i?i t 2, derived 
from the OPE and QCD perturbation theory, the following corrections have to be included: 
{%) electromagnetic radiative correction; {%%) effects of boost into the lab frame; {%%%) smearing 
due to detector momentum resolution. To take {%) into account, following the CLEO analysis, 
we used the resummed photon radiation corrections as given in Ref. [I8|. These corrections 



to R\ : 2 have very little sensitivity to subleading logarithms. To determine the corrections 
due to (ii), we assume that the B mesons are monoenergetic, with energy myus)/^ (the 
effect of the 4 MeV spread in the center of mass energy is negligible) . We found that the 
smearing due to the CLEO-II detector momentum resolution [TJjj , and the 50 MeV binning 
of the data has a negligible effect on R\^- 

Including the leading nonperturbative corrections of order (Aqcd/^) 2 [HH and the 



order a s corrections p0| , the theoretical expressions for Ri^ are 
i?i = 1.8059 - 0.309 ^- - 0.35 - 2.32 A- - 3.96 - ° 
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with Ri in GeV. We have defined the spin-averaged B meson mass, % = (m B + ?>m B *)/A. 
The terms in the last two parentheses in each of eqs. (|5a|) and ( |5b| ) represent the effect 
of the electromagnetic radiative correction, and the effect of the boost into the lab frame, 
respectively. The parts of these corrections proportional to A12 are negligible. Eqs. (0) 
correspond to electrons; for muons the electromagnetic correction is smaller. Corrections to 
eqs. (H) of higher order in a s and Aq CD /m6 will be discussed later. Even though there are no 
nonperturbative corrections to R\^ of order eqs. ([5]) contain terms proportional to 

A/m B . These arise since we reexpressed the heavy quark masses in terms of hadron masses, 
using m Q = m M - A + Ai/ (2m M ). 

To compare the above theoretical expressions with data, we need to discuss the ex- 
perimental uncertainties. We use the single tagged data to extract -R1.2, and correct for 
the effects of secondary leptons |2l| using the double tagged data. The central values are 
Ri = 1.7831 GeV and R2 = 0.6159 (the corrections from the secondaries are 0.0001 and 



0.0051 respectively), while the correlation matrix of the statistical errors [|T^] is 
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/3.8xl(T 6 6.0x10" 
V(R 1 ,R 2 )=[ 1 . (6) 

\6.0xl0- 6 1.7xl0~ 5 

The largest part of these uncertainties is due to the errors in the secondary lepton spectrum 
from the double tagged data. Estimating the systematic errors is more complicated. These 
uncertainties in the lepton spectrum can be divided into two classes: there are additive cor- 
rections, like backgrounds that are subtracted from the data; and there are multiplicative 
corrections, like those in efficiencies. The total systematic uncertainty in the CLEO mea- 
surement of the semileptonic B decay branching fraction is about 2%. However, only a small 
fraction of these uncertainties affect the shape of the lepton spectrum above 1.5 GeV ||22|| . 



In this region the uncertainties in the backgrounds are small, and the efficiencies have fairly 
flat momentum dependences. While the uncertainties in the electron identification and in 
the tracking efficiencies are the dominant sources of systematic error in the semileptonic B 
branching fraction, they are expected to affect at a much smaller level. We estimate 



that the systematic uncertainties in R 12 are not larger than the statistical errors pzfl , al- 
though a complete analysis of these can only be carried out by the CLEO Collaboration. 
For this reason, and since the statistical errors can be included into our analysis exactly, the 
experimental uncertainties we shall quote will be the statistical ones only. 

The comparison of the theoretical predictions in eqs. (|^) with the CLEO data is shown 
in Fig. 1. The steeper band is the constraint from R 2 , while the hatched one is that from Ri. 
The widths of the bands represent the la statistical errors, while the ellipse shows the la 
allowed region in {A, Ai}, after correlations between Ri and R 2 are taken into account. This 
region corresponds to A = 0.39 ± 0.11 GeV and Ai = -0.19 ± 0.10 GeV 2 . The la allowed 
region in Fig. 1 lies outside but not far from the region allowed by a recent analysis based 



on moments of the hadron spectrum |T2 



In Fig. 1 we set |\4&/Vd,| = 0.08. The extraction of this value is model dependent, and 
therefore has considerable uncertainty. If IV^/V^I = 0.1 then the center of the ellipse in 
Fig. 1 would move to A = 0.42 GeV and A x = -0.19 GeV 2 . We used a s = 0.22, corresponding 
to the subtraction scale mj. The sensitivity of our results to this choice of scale is weak; 




A [GeV] 

FIG. 1. Allowed regions in the A — Ai plane for R\ and i?2- The bands represent the la 
statistical errors, while the ellipse is the allowed region taking correlations into account. The star 
denotes where the order (AQco/w-b) 3 corrections discussed in the text would shift the center of the 
ellipse. 

changing a s to 0.35 moves the central values to A = 0.36 GeV and Ai = —0.18 GeV 2 . 

To plot Fig. 1 we used the data corresponding to electrons only, as we suspect that 
the systematic uncertainties in the (single tagged) muon data may be larger (for example, 
the muon detection efficiency is strongly energy dependent below 2 GeV). The latter data 
set, nevertheless, yields a consistent determination of A and Ai, giving central values A = 
0.43 GeV and Ai = —0.21 GeV 2 (to subtract secondaries we used the double tagged electron 
data). 

Theoretical uncertainties in our determination of A and Ai originate from the reliability 
of quark-hadron duality at the scales corresponding to the limits in the integrals defining 
Ri,2, from order (AqcdA^) 3 corrections, and from higher order perturbative corrections. 
Concerning duality, note that Ei > 1.5 GeV and 1.7 GeV (in the lab frame) correspond 
to summing over hadronic states X with masses below 3.6 GeV and 3.3 GeV, respectively. 
These scales are likely to be large enough to trust the OPE locally. This is supported by 
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the fact that a modified ratio that differs from i? 2 only in that the integration limit in the 
numerator is changed from 1.7 GeV to 1.8 GeV yields a parallel band that overlaps with 
that corresponding to _R 2 . Using this variable and R\, the central values for A and Ai 
become A = 0.47 GeV and Ai = —0.26 GeV 2 . (The assumption of local duality becomes 
less reliable using 1.8 GeV.) For higher moments |23| theoretical uncertainties increase, and 
they are sensitive to an almost identical combination of A and Ai as the first moment, 
R\. For example, the normalized second moment (with Ei > 1.5 GeV) gives a band that 
overlaps with that from Ri, and together with R2 yields the central values A = 0.39 GeV 
and Ai = —0.19 GeV 2 . Perturbative corrections of order a 2 and order a s (AQco/^fe) 2 have 
not been computed. The order a 2 corrections are likely to decrease the magnitudes of A 
and Ai. 

The recently calculated order (AQ CD /m fe ) 3 correction to the differential decay rate |24 



is parametrized by two matrix elements, p\ and p 2 , and time ordered products of local 
operators. These also modify the relation between the quark and hadron masses. Neglecting 
p 2 (that is expected to be small |Tj|) and the time ordered products, the ellipses in eq. (^) 
is replaced by pi/(2mg) 2 ; and A 2 is determined by m B , — m B = (2A 2 /7?7b)(1 + A/me). 

To estimate the possible size of these (AQ CD /m b ) 3 corrections, we calculated the effect 
of the terms proportional to A 3 , AAi j2 and p\ on the central values of A and Ai. The 
numerical value of p\ can be estimated using the vacuum saturation approximation, p\ = 
(27ra s /9) mg f B ~ (300 MeV) 3 (we agree with Jl3| and disagree with |25) on this result). 
Reexpanding i?i >2 to order (Aqcd/^) 3 , we get 

SR X = -(0.4A 3 + 5.7AAi + 6.8AA 2 + 7.7pi)/m| , (7a) 
SR 2 = -(1.5A 3 + 7.1AAi + 17.5AA 2 + l.Sp^/m^ . (7b) 

Solving the constraint imposed by the experimental data, our central values (the center of 
the ellipse in Fig. 1) are moved to the point denoted with a star in Fig. 1, A = 0.35 GeV 
and Ai = —0.15 GeV 2 . We think that the terms displayed in eqs. (0) are likely to be the 
dominant order (Aqcd/^6) 3 corrections. Uncertainties in the values of p 2 and in the time 



ordered products are expected to have a numerically smaller effect on R\^- Therefore, the 
(Aq C d/^&) 3 corrections will probably reduce the magnitudes of both A and Ai. 

A is not a physical quantity, and has a "renormalon ambiguity" of order Aqcd I^H • The 
perturbative expression for the semileptonic decay rate in terms of the b quark pole mass rrib 
is not Borel summable, and neither is the perturbative expansion of the MS mass m^m^) in 
terms of m&. These ambiguities cancel if A (or equivalently the b quark pole mass) extracted 
from the differential semileptonic decay rate is used to get the MS mass. Consequently 
one can arrive at a meaningful prediction for m&(m&). It is fine to introduce unphysical 
quantities like A, as long as one works consistently to a given order of QCD perturbation 
theory and the expansion in inverse powers of the heavy quark masses. Since the final 
results one considers always involve relations between physically measurable quantities, any 
"renormalon ambiguities" arising from the bad behavior of the QCD perturbation series at 



large orders cancel out p7 1 . 



In summary, using CLEO data on the lepton spectrum from inclusive semileptonic B — > 
X I Vi decay, we obtained the values of the heavy quark effective theory matrix elements, 
A = 0.39 ± 0.11 GeV and Ai = —0.19 ± 0.10 GeV 2 , where the uncertainty corresponds to the 
la statistical error. These imply at order a s that 

. . fBr(B -> XJPA 1.54ps\ 1/2 

W=°- M1 ( 0.105 r B ) ■ < 8) 

The largest uncertainty (~ 4%) in this determination of \V c b\ is due to higher order correc- 
tions in a s {e.g., the order a 2 s f3o correction to the inclusive semileptonic B decay rate | 28[ 
increases the factor 0.041 in eq. (§) to 0.042). 

The difference between the bottom and charm quark pole masses is free of renormalon 
ambiguities (at order Aq CD ). We find m b — m c = 3.37 ± 0.02 GeV, where the uncertainty is 
the la statistical error. The MS quark masses are related to the pole mass via mQ(mg) = 
mQ[l — Aa s {mo) / '(Sir) + ...]. Using this and our values for A and Ai, we obtain at order a s 

m b (m b ) = 4.45 GeV , m c (m c ) = 1.28 GeV . (9) 
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Order a 2 terms in the relation between ttiq and mg(mg) reduce these MS heavy quark 
masses by about 230 MeV and 180 MeV, respectively. Of course, it is not consistent to use 
these corrections as order a 2 terms have not been included into our determination of A and 
Ai through i? 12 . Recent determinations of the MS b quark mass using lattice QCD are 
Wi&(m&) = 4.17 ± 0.06 GeV and m^m?,) = 4.0 ± 0.1 GeV |[29|| . Given the uncertainties in 
eq. (|^) from order a 2 effects (and since the lattice results only contain the one-loop matching 
between lattice and continuum perturbation theory), our result for m^rrib) is consistent with 
these. 

The bands in Fig. 1 corresponding to R\^, are almost parallel. Hence, even small correc- 
tions to eqs. (H) can significantly affect our determination of A and Ai. It would be useful 
to have constraints that depend on very different combinations of A and Ai. The moments 
of the photon spectrum in inclusive B — * X s 7 decay fill can provide such information. 

Some improvements on the analysis in this paper are possible. The dominant (order a 2 /3o) 
part of the order a 2 corrections should be readily calculable f2~8fl . The order a s (Aq CD /mfe) 2 
corrections would set the scale for A2, and could also yield sizable corrections to the co- 
efficients of Ai in eqs. (^). A better knowledge of the order (AQ CD /m b ) 3 corrections and 
a model independent determination of \V u b\ would also reduce the uncertainties. The sys- 
tematic experimental errors on i? 12 should be evaluated by the CLEO Collaboration, and 
hopefully in the future both the systematic and the statistical errors can be reduced by a 
larger sample of double tagged semileptonic decays. Such a sample may also permit the use 
of variables analogous to -Ri^, but with smaller lower limits of integration, thereby reducing 
uncertainties associated with the validity of quark-hadron duality. 
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